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Abstract 

The Das-Mathur-Okubo (DM0) sum rule for the polarizabihty of charged 
pions is evaluated for the Nambu-Jona-Lasinio model Lagrangian in both its 
minimal and extended forms. A comparison is made with the results obtained 
using the same sum rule from chiral perturbation theory (xPT), approximate 
QCD sum rule calculations, explicit calculations on the lattice by Wilcox, and 
using the semi-empirical Kapusta-Shuryak spectral densities. The xPT results 
from Compton scattering are also given. We point to a delicate cancellation 
between the intrinsic and recoil contributions to a^± in the DM0 sum rule 
approach that can lead to calculated polarizabilities of either sign. 

PACS numbers: 11.30.Rd, ll.lO.Lm, ll.lO.St 
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The electric and magnetic polarizability coefficients and /^tt of pions are fundamental 
constants of strong interaction physics and their calculation is an important testing ground 
for QCD or effective models thereof p|-p!0|, since these coefficients can be measured exper- 
imentally. There is however, a reasonable amount of confusion as to their values, and a 
review of the current status of theory and experiment is given in []TT|. While a-,^ and (5-,^ can 
be identified directly from the soft photon limit of the pion Compton scattering amplitude, 
such calculations are extremely tedious, even in lowest order To obtain a physical 

picture and common understanding of the values obtained by several approaches, it is far 
simpler to use the Das-Mathur-Okubo (DM0) current algebra sum rule [O, and we exam- 



ine this here. We also accept the constraint Q;,r + /^tt = 0, valid for chiral pions |T3[, and 
concentrate on the electrical polarizability of charged pions. For charged pions, Holstein 
has shown that the DM0 sum rule offers an alternative route for calculating the intrinsic 
contribution to a.„± by recasting it asQ 



that together with the center of mass recoil contribution |T^ determines a.,^±: 



In these formulae, py^A are the vector and axial vector spectral densities, m^, (r^) and 
the pion mass, radius squared and decay constant, and a ~ 1/137 is the fine structure 
constant. The summation in the defining expression for a*"*'" in Eq. (|I|) runs over all electric 
dipole excitations |i) of energy Ei that are connected to the pionic ground state by the 
dipole operator d, and the minus sign on the second term anticipates the fact that a*"**" < 
for pions. Physically, this comes about from an interplay between negative energy 
intermediate states entering the sum in Eq. (|1|), and the contributions of disconnected photon 



and pion decay diagrams that have to be subtracted out ||16 . 

The DM0 sum rule result for a*"^*'' is valid to O(p^) in the language of chiral perturbation 
theory (xPT) [jl^ due to the approximations made in its derivation. This can be seen as 
follows. By inserting the definitions of the Gasser-Leutwyler coupling constants /s, of xPT 

1 - /""^ ds 1 - 

<r^>=Y^(/6-i); I ^[pW^)-pa(3)] = ^(/5-i) (3) 

into Eqs. (|lD and (D, one retrieves the xPT expression at one loop level (i.e. also to (9(p^)) 
for as has also been extracted from Compton scattering 

OL - - 



^The integral on the right carries an extra 1/2 in |14| due to the conventional QCD lattice 
normalization of the spectral densities used there. 
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This takes the numerical value 2.7 x 10~^ fm^ using experimentally extracted values of the 
li, and is a direct prediction of QCD. We may thus only make a comparison with xPT to 
this order in using the sum rule. Note that a higher order calculation of has been made 
via Compton scattering 0. This calculation involves the computation of over 100 Feynman 
diagrams, and lowers the final result from = 2.7 to 2.4 xlO~^ fm'^. 

The DM0 form for a*"**^ in Eq. (0) has the particular advantage that this contribution 
to a^ri can now be estimated using QCD sum rules as in or calculated directly on the 
lattice as in [Q. Here, we also implement the DM0 calculation of a*"**^, as well as the 
recoil contribution to a^ri, using a model chiral Lagrangian given below to describe the 
quark dynamics, and we then compare this with QCD sum rule results, calculations on the 
lattice, chiral perturbation theory, and also a result obtained using the empirical Kapusta- 
Shuryak [|I8| spectral densities to evaluate the DM0 sum rule. Recent data from the ALEPH 
collaboration [l^j2^ for the spectral densities is also shown. We base our discussions on the 
extended Nambu-Jona-Lasinio (ENJL) Lagrangian [^,0 which explicitly includes vector 
and axial vector degrees of freedom: 



-ENJL 



(5) 



Here ip = {u, d) is the SU (2) quark flavor doublet, r Pauli isospin matrices, Gi, G2 coupling 
constants, and m a common current quark mass that explicitly breaks chiral symmetry. 
Together with a regulating cutoff A, these parameters are fixed by reproducing the pion 
mass and decay constant, quark condensate density, and the axial coupling constant of the 
constituent quarks. 

The spectral densities pj{s) = — (l/47r)/mHj(s) at four-momentum transfer squared 
s = q"^ that enter into Eq. (|I]) are obtained from the two point polarization functions of the 
vector and axial vector currents J^{x) = V"^(x) or ^^(a;) of isospin index a and have the 
tensor structure H^j,.„j(s) = Ilj{s){g^u — q^^qv / q^)5ah- These amplitudes are purely spin one 
("transverse") in the chiral limit due to current conservation, d'^J^{x) = 0. Their difference 
satisfies the unsubtracted dispersion relation p3|j2^ 



dt 



t 



(6) 



We now specialize to the case where the dynamics of the are assumed to be governed by 

Cenjl- Working to leading order in A^""^ P3|P^ , both the lowest order one- loop expressions 
Hv,Ai ^"^^ resummed expressions for the polarization functions Hy,^ as given by their 
Bethe-Salpeter equation, can be explicitly shown to satisfy the unsubtracted dispersion 
relation, Eq. (|^), under Pauli- Villars regular iz at ion0 Using this dispersion relation, the 
required integral in Eq. (|I]) can be rewritten as 



dt 



t 



(7) 



^In which case the associated spectral density difference in Eq. ^ mimics |27] the QCD asymptotic 
behavior pv{t) - pA{t) ~ C(s"^) 
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The difference of polarization amplitudes in the ENJL model is found by direct calculation 
to be given by the relation (Ward identity |28|) 



^[flAis) - Ily{s)] = f^Fp{s)Fv{s)FA{s) (8) 

in terms of the pion electromagnetic form factor Fp{s) of the minimal NJL model {G2 = 0) 
and the vector and axial vector form factors Fv{s) and Ga^s) of the constituent quark 
currents. Closed expressions for the latter two form factors are obtained by summing the 
associated Bethe-Salpeter equation to all orders in the coupling 2G2- One finds 

Fv{s) = [1 + 2G2Uv{s)]-'; Ga{s) = [1 + 2G2Ua{s)]-' = gAFA{s) (9) 

in terms of the irreducible one-loop polarization functions Ylv,A- Here qa = Ga(0), or I—qa = 
8G2f^, is the axial coupling constant pT H33|j35[] . The derivative on the right hand side of 
Eq. (1^) that is required to evaluate the integral explicitly is determined by the behavior of 
the form factors at vanishingly small momentum transfers. This can be found in terms of the 
p and ai meson masses mv,A and their common coupling constant Qp to the quarks that can 
be identified from the poles and residues in the corresponding resummed propagators for the 
vector and axial vector modes —iDv,A{s) = {—2iG2Fv{s),—2iG2GA{s)} after introducing 
low-energy expansions for the irreducible polarizations that they contain: inv'(s) ~ 
—Qp'^s and |n^(s) ~ ''Qp^is — 6m^). Here m is the self-consistent quark mass of the ENJL 
model. One finds m\ = Qm? + my = rriy/gA and my = (7^/86^2 = dpfn/i^ ~ 9a)- The 
latter result follows by eliminating G2 in terms of its expression for qa given above. This 
last form for my reduces to the KSRF relation for the choice qa = 1/2. All of these 



low-energy relations also follow immediately by working with the bosonized version of the 
ENJL Lagrangian from the start [Q. Introducing the same expansions of the irreducible 
polarizations into the denominators of the form factors defined above in Eq.(P), one obtains 
Fy(g2) = 1 + q^/ml + ■ ■ ■ and Fa^Q^) = 1 + gV"^A - 9a{1 - gA)qV (Sn^ f^) + ■■■ after also 
recalling that Fp{q^) = 1 -|- fi'AQ'^/(87r^/^) Using this information the value of the 



integral in Eq. (|^ becomes 

r '^[Pvis) - Pa{s)] = |\ + 4 + 4, (10) 
JO ■s"' my m\ 

The pion radius parameter is found from the expansion of the ENJL pion form factor 
FM^) ^ Fp{q^) to be 

<^'>=7% + 4 + 4 (11) 
47r2/2 ml m\ 

for Nc = 3. Referring back to Eqs. ([l|) and (@), one sees that an exact cancellation is going to 
occur between the meson mass terms in the recoil and intrinsic contributions to the electric 
polarizability. The final result for the ENJL polarizability as calculated using the DM0 sum 
rule thus reads 

^9a 
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To the order 0{N^) and 0{p^) we are working,^ this result depends only on the axial 
coupling constant of the constituent quarks besides the pion parameters. The lesson from 
this is that there is a delicate balance between the recoil and intrinsic contributions, so 
that at a more fundamental level, one should probably compute both terms to the same 
level of approximation, instead of relying, for example, on the known experimental pion 
radius squared for the first term. As in this can give negative values for a.,^±. If, 
in addition, we insist that the KSRF relation P3 should also hold between the meson 
parameters generated by the ENJL model in its low energy regime, then we have seen thatQ 
= 1/2, and Eq. (|T^ becomes a parameter-free result for the charged pion polarizability, 
of magnitude a.,^± ^ 1.5 x 10~^ fm^ for the physical constants. This value is listed in Table 
I. 

The expression for as obtained from Compton scattering calculations using the NJL 
model |0,|lO| also reduce to Eq. (|l^) (with 5'a=1) when evaluated to leading order OijP'). 
This confirms that the integrals over the model spectral densities behave correctly and is 
an important self-consistency check, since these spectral functions give the density of states 
in the qq continuum of the non-confining Cenjl-, and not that of the physical p ^ titt and 
ai —>■ TTp continua. Notwithstanding this unphysical feature, the sum rules continue to be 
obeyed. The reason for this is that the model polarization amplitudes satisfy unsubtracted 
dispersion relations like Eq. and the Ward identity of Eq. (|) under gauge-invariant 
regularization schemes (such as Pauli-Villars or dimensional regularization). This in turn 
allows one to re-express the sum rules in terms of the zero momemtum transfer values of 
the model polarization amplitudes and their derivatives that are not sensitive to the non- 
confining nature of the model. Fig. |I] compares the behavior of the vector minus axial 
vector spectral density difference for the ENJL model with (a) the semi-empirical fits of 
Kapusta and Shuryak [0, and (b) recent data from the ALEPH collaboration []T^^. Both 
theoretical sets of curves satisfy Weinberg's first sum rule |^3[. The ENJL model density 
difference accomplishes this by having considerably smaller and wider resonant peaks at the 
(model) vector masses my = 793 MeV and rriA = 1122 MeV. Note, however, that the value 
that we obtain for a^ri using the Kapusta-Shuryak densities is negative as is seen in Table I. 
One sees that the vector peak of the ALEPH data is well-modelled by the phenomenological 



form of Ref. [|18[. At large s, the uncertainties in the data are too large to make it possible 
to evaluate the sum rule precisely, for s — > oo [pO| . 

We next return to the NJL model and examine the sub-leading in contributions to 
Eqs. (|TI|) and (|TU|) and show that they are not small. We illustrate these corrections in the 
case of the minimal NJL model {G2 = therefore g^ ^ 1 and my,m\ — » 00) that only 
supports 71 and a meson modes. Such a calculation has already been performed for the pion 
radius to this order, and one has found ||39|j40[| 



^Since G2 ~ A^^ ^ in order to preserve the correct scaling properties of quark and meson masses, 
and ~ A^c; one sees that 1 — gA = 8^2/^ ~ N^!. The same result follows from counting rules 

^ Estimates of gA using the Adler-Weisberger sum rule |37,38|, suggest a much larger (~ 0.8 — 0.9) 
value for gA- 
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which is necessary for evaluating the recoil contribution to q;^± . Note that the /^^ which en- 
ters this expression (via a Goldberger-Treiman relation that is itself correct to 0{1/Nc)) 
refers to a pion weak decay constant that is correct to 0{1/Nc). This comes about because 
(r^) follows from differentiating an NJL electromagnetic form factor for the pion that in- 
cludes 0{1/Nc) corrections, and therefore contains a nqq coupling constant g-^qq that is also 
correct to this order. Otherwise the modified form factor does not normalise properly to 
unity at g2 = 0. This is demonstrated explicitly in P0| . 

A similar analysis can be made for the spectral density function integral. This is briefly 
sketched here. Corrections arise from two types of diagrams : (i) single meson loop dress- 
ings of the irreducible one-loop quark polarization diagram, and (ii) meson pairs exchanged 
between p and Oi vertices. While type (i) is essential for maintaining the transverse nature 
of the vector polarization diagrams to 0{N^), the subset (ii) of irreducible two meson loop 
diagrams is singled out by the fact that their derivatives at = contain chiral logarithms 
Pl| that diverge as —>■ due to the pion pole in the propagator of the pion. So 
they dominate over the derivatives^ of type (i) in the chiral limit. Denoting the meson pair 
exchange diagrams with the coupling to the external channel factored out by —ill'^^{q'^) 
and —iIl'J^]J"^{q'^) respectively, one has 

i-l9,?[]Ti';.Hq')] = J ^r^^^(/,/ + g)^.((^ + g)')^.(/')rr(^ + ?,0 (14) 

for the p channel. A similar expression holds for the ai channel. The leading contribution 
in chiral logarithms to this is found most directly by using the bosonized version p5| , pl| of 



the Lagrangian (^) to identify the effective meson vertices T'^'^ and T^^g'^'^. These vertices 
are generated by the Lagrangians 

SCp^^ = -gpeabcT^ad^TibPc \ ^l^a^-Ka = gpCrC^Sl ■ O^lf + TT ■ 9^5'') (15) 

where (a, vr, p, ai) are isoscalar and isovector meson flelds and Qp is the common vector meson 
coupling constant. Filling in the details, one flnds the convergent contribution 



- nF-(.)].,o ^ -L /' dead - a) In { '"''^ " °> + '"'° } 

AStt"^ ^ mi 



to be added to Eq. (|^), after projecting out onto the transverse parts of the meson loop 
polarization diagrams. Here and nia = {rri^ + Arn?)"^ are the tt and a meson masses 



^In fact the type (i) diagrams are strictly constants i.e. are independent of the external meson 
momentum in the bosonized version of the NJL model and thus give no contribution at all to the 
difference of polarization derivatives in Eq. (^) in that case. 
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of the bosonized NJL model [0. The revised value of the integral ([T0|) for the NJL case 
{gA = 1) now reads 

I ^[p^(s)-p^(s)](quark+mesonloops) = ^-^(ln^ + ^)+0(m^). (17) 

In Eqs. ( P!7[ ) and ( ]TB| ) we have also retained chiral 0{N^) corrections, which are pure numbers 
in addition to the chiral logarithms that arise from the meson loop, in order to be completely 
consistent in Nc power counting. From Eqs. (^, one obtains values for the two scale- 
independent Gasser-Leutwyler coupling constants I5 + ln(m^/yU^) = 7 — 5/6 = 37/6 and 
lQ + hi{ml/ fj?) = 13 — 17/6 = 61/6, (here evaluated at the natural scale /i = rria ~ 500 — 600 
MeV of the NJL model) that are in the form given by chiral perturbation theory [T7 . 



The chiral logarithms cancel in the revised value of their difference which now becomes 
^6 ~ ^5 = 4 in Eq. (^), down from 13 — 7 = 6, due to the remaining 0{N^) mass-independent 
corrections that the Us now contain. This is to be compared with the empirical difference 
of 16.5 — 13.9 = 2.6. The 0{N~^) meson loop corrections are thus large and negative in this 
case, leading to a considerable decrease in the predicted NJL value for Table I. 

To summarize, the charged pion polarizability is examined in the NJL and ENJL models 
and compared with the results obtained in xPT, which serves as a benchmark of QCD. 
For completeness, the Compton scattering result for this quantity at two loop order is also 
quoted, hile this does not suffer from any obvious cancellation problems, it involves the 
computation of more than 100 graphs. The NJL model can be directly be compared with 
xPT to one loop order, in that the formal expression (4) holds in the NJL model with low 
energy constants evaluated within the model. This, in turn, leads to a much higher value of 

than does xPT, unless one includes the next order meson loops in the NJL model 
expansion. In the ENJL model, to lowest order, the final result is dependent on the value of 
gA- In an analysis according to the division into intrinsic and recoil contributions, one notes 
that the intrinsic contribution to the polarizability changes by about a factor 3 from the 
minimal NJL calculation {G2 = 0) to using the empirical densities of Kapusta and Shuryak. 
The negative results for a-,r± in the latter well as for the lattice calculations, reflect 

the delicate balance between the two relatively large recoil and intrinsic contributions to 
Q;,r± that are of opposite sign. This emphasizes the need for calculating both terms in the 
sum rule approach to the same level of approximation. The ENJL result (—11.0) for a™*'' 
is itself nicely bracketed by the QCD sum rule result (—9.6) and chiral perturbation theory 
(—12.6), all in units of 10~^ fm^, having been considerably increased over the minimal NJL 
case (—5.9) by the presence of the vector and axial vector modes. We have set gA = 1/2 
in the ENJL case in order to comply with the KSRF relation. However, even allowing gA 
to fall through its full range 1 > gA ^ (while scaling m like gA^^'^ to keep the same /^) 
gives values of the integral that only range between 12.7 and 27.5. In the latter instance the 
intrinsic and recoil contributions would simply cancel, leaving a polarizability coefficient of 
zero. 
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FIGURES 



FIG. 1. The vector minus axial vector spectral densities as calculated from the ENJL model 
(solid curve) , and from the semi-empirical Kapusta-Shuryak expressions of Ref . |jl^ (broken curve) 
versus momentum transfer squared. The dotted points with error bars are recent experimental 
values taken from the ALEPH collaboration. The Pauli-Villars scheme has been used to regulate 
the ENJL polarization amplitudes. This introduces two further artificial thresholds (only the first 
one is shown) at equally spaced intervals 4A^ beyond the unphysical qq threshold that lies at 
4m2 = 0.28 GeV^ for the particular ENJL parameter choice d = 2.47 GeV-^, G2 = 3.61 GeV-2, 
and a regulating cutoff of A = 1.06 GeV, used in this figure. 
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TABLES 



TABLE L Breakdown of intrinsic and finite size contributions to the polarizability of charged 
pions in units of 10~^ fm'^ (1 fm ~ 197.2 MeV~^) as obtained from various theoretical approaches, 
together with the available experimental data. The original Serpukov II and Mark II data analyses 
quote (a^± ± (3^±) without assuming the constraint a^± + /3^± = 0. We have simply averaged the 
statistical and systematic error bounds when extracting a^±. Input: m,r = 139.6 MeV, f-,^ = 93 
MeV, and a typical quark mass m = 324 MeV from solving the gap equation |21] for the ENJL 
case, giving my = 793 MeV and = 1122 MeV; empirical values = 13.9 ± 1.3, = 16.5 ± 1.1 
for the xPT case |l^, and the measured value < >= 0.439 fm^ |4^ except where derived from 
a model. 
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NJL [Eqs. (|10|,|1D, gA = l] 

NJL+0(iV-i) [Eqs. (0,0)] 

QCD sum rule g 

ENJL [Eqs. qa = 0.5 (from KSRF)] 

Lattice calc 







Semi-empirical densities [18] 
Lebedev 



Serpukov I [44| 



Serpukov II [45| 



Mark II |6| 

77 — !■ vrvr measurement 



20 ±12 
6.8 ± 1.4 
8.5 ±4.8 ±3.5 
2.5 ±0.5 ±0.02 
2.2 ± 1.6 
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